The in-medium scale evolution in jet modification 
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The in-medium modification of the scale dependence of the fragmentation function in dense mat- 
ter, brought about by higher twist corrections to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations, is derived. A phenomenologically motivated resummation is out- 
lined which incorporates the next-to-leading twist single gluon emission kernel along with the vac- 
uum emission kernel and provides an in-medium virtuality evolution of the final fragmentation 
function of a hard jet propagating through dense matter. The concept of a fragmentation function 
is generalized to include a dependence on distance travelled in the medium. Following this, numer- 
ical implementations are carried out and compared to experimental results on the single inclusive 
suppression observed in Deep-Inelastic scattering (DIS) off a large nucleus. 
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With the commissioning of the Large Hadron Collider 
(LHC) heavy-ion experiments, the study of hard probes 
in deconfined matter will have entered the deep pertur- 
bative domain: one expects jets with energies of sev- 
eral hundred GeV to be produced, plough through the 
dense deconfined matter and fragment far outside the 
medium [l|. With the start of the future Electron-Ion 
Collider (EIC), a similar step will be taken in the study of 
cold nuclear matter with hard jets. Current experiments 
at the Relativistic Heavy-Ion Collider (RHIC) @, [1] and 
at DESY (HERMES) [j| have just begun to probe the 
perturbative regime. At RHIC this is believed to start 
above a rather high transverse momenta pt ~ 6 — 7 GeV 
(due to the expanding medium) In DIS, the extent of 
the perturbative region depends on the energy (y) lost by 
the electron, the fraction carried by the detected hadron 
and the size of the struck nucleus [6| . 

In the case of single-hadron-inclusive suppression in 
DIS on a large nucleus A, the application of perturbative 
methods is based on the factorized cross section to pro- 
duce a hadron with a fraction z of the photon light-cone 
momentum q~ Q: 



da h (Q 2 
dz 



(1) 



In this equation, Fi{xi) is the parton distribution func- 
tion for a parton % with momentum fraction x%, cr^j is 
the hard cross section with the virtual photon to produce 
parton j and D is the medium modified fragmentation 
function to fragment into the hadron h after undergoing 
multiple scattering in the nuclear medium 0). 

The multiple scattering of the hard parton in the 
medium gcnerically has two parts, a perturbatively cal- 
culated contribution which represents the scattering and 
gluon radiation from the hard parton, and the non- 
perturbative distribution of the soft gluons off which the 
hard partons will scatter . Both these parts in combi- 
nation with the standard vacuum shower of gluons lead- 



ing to the non-perturbative fragmentation into observ- 
able hadrons are included within the definition of the 
medium modified fragmentation function. The three dif- 
ferent factorized parts in Eq. (JTJ) will necessarily depend 
on the factorization scales (//j for the initial state and \i / 
for the final) which in Eq. ([I]) have been chosen as the 
hard scale Q 2 . So defined, the medium modified frag- 
mentation functions will reduce to the standard vacuum 
fragmentation functions D(z,Q 2 ), as the extent of the 
medium in the final state is reduced. 

As the energies of the process are increased, all three 
factorized functions in Eq. Jl} change with Q 2 . While 
volumes of work exist on the perturbatively calculable 
scale dependence of both the hard part and the nuclear 
structure functions (as well as vacuum structure and frag- 
mentation functions) 13, II , IH , little attention has been 



paid to the scale dependence of the medium modification 
of the fragmentation functions. In the case of DIS on 
a nucleon or in p-p collisions, i.e., in the absence of an 
extended medium, the vacuum fragmentation functions 
have a well known and perturbatively calculable depen- 
dence on the scale Q 2 of the process given by the DGLAP 
evolution equations (Toj : 



dD h q {z,Q 2 ) _ a ,(Q2) r^dy 
d\og(Q 2 ) 2tt J z y 



-P q ^i(y)DN-,Q 2 ).(2) 
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In the above equation, a is the strong coupling constant 
and P q ^i represents the probability for a parton q to 
split into an parton i which carries a fraction y of the 
original momentum and finally fragments into the hadron 
h (a sum over i is implied). To date, the analogue of 
such an evolution equation in a medium has not been 
given. It is the object of this Letter to carry out such 
an extension and to present the first numerical results 
of such an evolution in an extended dense medium. An 
alternative approach to this problem leading to a Monte- 
Carlo implementation has been proposed in Ref. 

It should be pointed out that the factorized form in 



2 



Eq. dp i s an assumption that has not yet been proven to 
hold [14[. So far such a factorized form has been shown 
to hold (at leading twist) in the absence of a medium. 
Factorization of the hard part from the initial distribu- 
tion has also been demonstrated to hold at all twist in 
the case of the totally inclusive cross section in DIS and 
up to next-to-leading twist for the single-hadron-inclusive 
cross section [3] ■ While it may reduce the rigor of the 
present approach, we persist with this form (in the ab- 
sence of a more accurate formalism), seeking justification 
in phenomcnological applications. While the LHC and 
EIC will engender energy scales that are higher by an 
order of magnitude, current experiments at RHIC and 
HERMES already sample a wide enough range of scales 
to provide a sufficient test for this formalism. 

Imagine that a quark jet is produced in a hard colli- 
sion with a large virtuality M 2 . This virtuality is then 
reduced by a sequence of partonic emissions. The emit- 
ted partons will possess a lower virtuality than the par- 
ent. These will then decay further into lower virtuality 
branches and the process will continue. The subsequent 
branchings and the development of the partonic shower 
may be calculated using pQCD as long as the virtual- 
ity at any vertex is large compared to Aq CD . Once the 
virtuality drops to some predetermined value fi 2 (a few 
times Aq C £)), a non-pert urbative object such as a frag- 
mentation function [l5| needs to be introduced. 

In order to incorporate the effect of a medium on this 
cascade process, the fragmentation function has to be 
generalized to include position dependence, i.e., 



D(z, M 2 



D(z,M 2 ,q-)\f_ 



(3) 



denotes the location of the production of a hard quark 
with large negative light cone momentum q~ and virtual- 
ity M 2 in an extended dense medium. The leading parton 
from the shower ensuing from such a quark is assumed 
to exit the medium at location (J. We will insist that 
the final exiting parton which fragments to produce the 
detected hadron possesses virtuality to 2 (^ ^qcd) that 
is much larger than the transverse momentum gained by 
the jet in traversing the medium (we chose this to be > 1 
GeV for jets punching through large nuclei). 

The virtuality drop from M 2 to m 2 may be achieved 
through any number of emissions. In the case of single 
emission without scattering, there is no dependence on 
location and the change in D{z, M 2 ) is given as 



M 2 1 



(4) 



The splitting function P q ^i(y) above contains the prob- 
ability for a quark q to radiate a gluon with momentum 
fraction y and transverse momentum l± (Io| . Virtual cor- 
rections which conserve unitarity arc implicitly included. 



The case of single gluon radiation accompanied by one 
scattering is given by the diagrams in Fig. [I] Only the rel- 
evant amplitudes are drawn. The gluons attached to the 
dark circles are meant to indicate gluons exchanged with 
the medium when the hard partons scatter off it. The 
light red (shaded) ellipses on the Feynman diagrams de- 
note the propagators that are considerably more virtual 
than the remaining parton lines. To the five diagrams 
included, one must add contributions with no scattering 
as well as contributions denoting the interference of dou- 
ble scattering with no scattering diagrams and virtual 
diagrams (denoted as v.c). The two diagrams in the 
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FIG. 1: Some of the diagrams included in the quark in- 
medium splitting function. See text for details. 



top row of Fig. [TJ where the initial parton with a large 
virtuality splits into two partons with lower virtuality, 
followed by a single soft scattering, is similar to the case 
in vacuum. The single soft scattering merely introduces 
a mild change in direction of propagation of the struck 
line. These diagrams are referred to as the vacuum-like 
contribution to in-medium induced radiation and their 
interference with the remaining diagrams as the vacuum- 
medium-induced interference contribution. The three di- 
agrams in the bottom row represent the case where the 
scattering with a hard gluon in the medium leads to a rise 
in the virtuality, immediately lost by gluon emission 0). 

The result of the interference of all such dia- 
grams [which leads to the Landau-Pomeranchuck-Migdal 
(LPM) effect] yields the "medium dependent" correction 
to the in-medium modified fragmentation function. If 
the drop in virtuality from M to m were to occur solely 
through single scattering and single emission, then @, 0], 



*J%\z t M* tq -)\% = %j 
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(5) 



In the equation above, p g {C) is the gluon density at (. 
The (— ) superscript on the locations have been dropped 
for simplicity. The tilde on the splitting function in- 
dicates that the color factor and virtual correction are 
somewhat different from the case of the vacuum (requir- 
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ing the final outgoing line, after a scattering and an emis- 
sion to be nearly on-shell restricts the range of y, see 
Ref. [H for details). In Eq. ([5]), £ represents the location 
of the scattering vertex of the hard parton off the medium 
in Fig. [1] Due to interference effects, the medium mod- 
ified fragmentation function is now a function of the jet 
energy q~ ; it is no longer universal and strongly depends 
on the details of the medium as introduced by the density 
factor p g (C) an d its space-time distribution. However, 
the functional dependence of the fragmentation function 
on the medium density is still universal and thus may 
be used to compute the medium modified fragmentation 
function in any medium where p g (C) is known. 

Given that each radiation leads to a drop in the vir- 
tually of the propagating parton, successive radiations 
are assumed to be strongly ordered in transverse momen- 
tum. While this is well known for the case of radiations 
in the vacuum, in a medium, the strong ordering of the 
virtualities will be broken if the parton encounters a hard 
scattering which will lead to a large transverse momen- 
tum radiation. At next-to- leading twist, as is the case 
in this calculation, this hard scattering can at most oc- 
cur at the location of the first rescattering (in the case 
of higher twist contributions as in Ref. [l6|, this hard 
scattering may occur at any later rescattering). In such 
an ordered scenario, we may write down the expression 
for the parton to lose virtuality from M 2 to m 2 in two 
ordered emissions as (we suppress q~ in the argument) 

M 2 1 Cf 

m 2 z 



2 - 2 cos 



ii(c-co 

2q-y{l-y) 
i . Cf 
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(6) 



Due to the drop in virtuality after the first emission, the 
location of the first scattering in the medium (whether 
soft or hard), is considered as the origin of the (\ inte- 
gration and thus of the interference pattern connected 
with the second scattering. Alternatively stated, this 
means that the above expression focuses only on ladder 
diagrams both in the scattering and in the emitted gluon 
sector. The insistence on this is simply based on the 
dominance of terms which contain a strong ordering of 
virtualities in successive emissions. 

In the strongly ordered or planar emission limit, the 
emission points for successive lower l± emissions are also 



strongly ordered. The extension to multiple emissions, 
may now be written down similar to Eq. ((6|). Differ- 
entiating the equation for D(z, M 2 )\^ f (containing mul- 
tiple emissions) with respect to log(Af 2 ) leads to the 
"medium dependent part" of the evolution equations for 
the medium modified fragmentation function, 

dD*{z,M*g-)\% a s fdy ? - 



aio g (M 2 



2tt J y 



n\[- y ,M*q 



The single scattering kernel K is given as 



K 



2ira s p(Q 

Nr. 



2 - 2 cos 



M 2 (C-Ci) 
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(7) 



(8) 



The full evolution equation for the medium modified frag- 
mentation functions will include contributions from pure 
vacuum splitting functions as well as contributions from 
gluon fragmentation functions which have a similar in- 
medium evolution. 

The simplest application of the formalism developed 
above is to compute the medium modified fragmentation 
function in the case of DIS on a large nucleus, where, at 
least, one hadron with a large forward momentum is de- 
tected in the final state. Experiments present the ratio 
of this with the vacuum fragmentation function at the 
same z and Q 2 , called the nuclear attenuation factor 0]. 
Experimental results for three nuclei (Ne, Kr and Xe) 
are presented in Fig. [21 Assuming single scattering and 
single in-medium emission, the medium modified frag- 
mentation function may be calculated using Eq. ([5]); the 
corresponding attenuation factor is the red dashed line 
in Fig. [21 In the case of multiple emissions, one will need 
to use the in-medium evolution equations of Eq. ([7} . In 
either case, the basic kernel K q - M 2(y, £) is identical. To 
calculate this kernel, we introduce a nucleon density in 
the large nucleus as an input Ansatz. Due to the simplic- 
ity in analytic calculations, we use a hard sphere density 
distribution (with radius Ra)'- p(C) = Po@{Ra — ICD- 
While Ra depends on the nucleus, po is a fit parameter 
dialed to obtain the best overall fit. We also approximate 
D(z)\^ s ~ D(z)\ J, which greatly speeds calculation. 

As may be seen from the comparison with the data in 
Fig. [21 the fit, for the case of single scattering and single 
emission (red dashed lines), with smaller nuclei such as 
A^e is adequate. The comparison, progressively worsens 
as one proceeds to larger nuclei. This is clearly seen in 
the case for Xe where there seems to be almost a differ- 
ent slope with z between the calculations and the exper- 
imental results. In some ways, this is to be expected; as 
one proceeds to larger nuclei, the possibility of multiple 
scattering and multiple emission increases and the results 
of a formalism which only included single scattering and 
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FIG. 2: Comparison of the experimental results on the nu- 
clear attenuation observed in the DIS on Ne, Kr and Xe 
with three different calculations of the ratio of the medium 
modified to the vacuum fragmentation functions. See text for 
details. The data points (and the bins of v and Q 2 ) are from 
Ref. [H. 



single emission in medium should show a systematic de- 
parture with nuclear size. This is to be contrasted with 
the medium modified fragmentation functions calculated 
using Eq. (|7|) (solid blue lines). One immediately sees a 
marked improvement in the comparison between the cal- 
culations and the experimental data. This represents the 
principal result of this Letter. 

While the inclusion of multiple emissions, each with a 
single scattering kernel leads to a marked improvement, 
there still is a small difference between the data points 
and the theoretical curves for the case of DIS on Xe at 
small z. One way to improve this is to include multiple 
scatterings per emission as was done for photon emis- 
sion in Ref. [l6j|. Such an effect is expected to become 
important for larger A, but will have to be left for a fu- 
ture effort. Yet another effect is the presence of a certain 
amount of hadronic energy loss [lH, in addition to the 
partonic energy loss calculated in this Letter. This ef- 
fect should be most dominant at lower z. To estimate its 
importance, we introduce an phenomcnological extension 
as indicated by the green dot-dashed line in Fig. [5] In 
this case the contribution from the gluon fragmentation 
functions is removed, i.e., hadrons emanating from the 
fragmentation of soft gluons radiated by the jet are com- 
pletely stopped in the medium. This latter calculation 
has the best agreement with the largest nucleus (Xe), 
where the agreement with data is remarkable. The ac- 
tual attenuation will lie in between the green dot-dashed 
and the blue solid lines. In each of the calculations, the 
gluon density was dialed to provide the best overall fit to 
all three data sets. This leads to a small spread in the 
medium gluon densities (po) required, which corresponds 
to a spread in the quark jet quenching parameter go (see 
Ref. [igj]) between 0.06 to 0.18 GeV 2 /frn. 



Applications of the methods presented in this Letter 
to the case of high transverse momentum hadron produc- 
tion in heavy-ion collisions has recently been presented 
in Ref. Q, where Eq. ([7]) was evaluated in a three di- 
mensional hydro-dynamically expanding medium. Due 
to the larger error bars in the data at RHIC, the im- 
provement of the present formalism over the medium 
modified fragmentation functions calculated in the single 
emission limit (in an identical medium, using Eq.© (l9() 
are somewhat less evident. 
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